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1 Introduction 



The theory of one-parameter semigroups provides a good entry into the study of 
, the properties of non-self-adjoint operators and of the evolution equations associ- 

^ ' ated with them. There are many situations in which such an operator A arises by 

■ linearizing some non-linear evolution equation around a stationary point. The sta- 

bility of the stationary point implies that every eigenvalue of the semigroup Tt = e"^* 
O '. has negative real part, but the converse is not true. This was vividly demonstrated 

^ I in a famous example of Zabczyk, in which the semigroup norm grows exponen- 

tially, although every eigenvalue of the operator in question is purely imaginary, 
jSl Th. 2.17], Pl]. One of the main points of this paper is to emphasize that similar 
phenomena occur for the so-called Schrodinger semigroups, which have extensive 
applications in quantum theory and stochastic processes. We will see that the long 
time behaviour of the norms of diffusion semigroups with self- adjoint generators 



^ ' may be entirely different for the and norms, although the generator has the 

c3 ■ same spectrum in the two spaces. In other words, growth bounds proved using 

the spectral theorem for self-adjoint operators may not generalize to the 'same' 
evolution equation acting in other Banach spaces, even when the other norm is 
physically more relevant than the Hilbert space norm. 

We mention in passing that in hydrodynamics Trefethen and others have estab- 
lished that pseudospectral methods may provide stability information unavailable 
by the use of spectral theory alone; see |2Z1I2H1E1- On the other hand Renardy has 
shown that in a number of hydrodynamic problems spectral analysis does indeed 
suffice to determine stability, [121 12011111 • 

There is an enormous literature studying the asymptotic behaviour of one-parameter 
semigroups as t — oo, |TT1, but as far as stability is concerned short time bounds 
on the semigroup norm are often more relevant: if ft = Ttf grows rapidly for some 
time, before eventually decaying exponentially, then the linear approximation may 
become inappropriate before this decay comes into effect (or would do under the 
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linear approximation). The fact that the short time and long time behaviour of a 
semigroup may be quite different is physically very clear for the convection-diffusion 
operator on a bounded interval or region, [iH], |HD1 p. 16-19]; see also 0- In this 
case the underlying cause is the non-self-adjointness of the operators concerned, 
which act in a Hilbert space. 

Our goal is to obtain information about the short time behaviour of the semigroup 
from norm bounds on the resolvent operators - closely related to the pseudospectra, 
for which efficient computations are now available, EZl Ell EUl IS2] • We succeed 
in obtaining lower bounds, not on the semigroup norms themselves, but on certain 
regularizations, defined in the next section. We also show (Theorem that it is 
not possible to obtain similar upper bounds from numerical information about the 
resolvent norms, however accurate this information may be. Both of these facts 
are completely invisible if one only looks at the spectrum of the relevant operator, 
which is of limited use for stability analysis. 

Some of the results in this paper are already familiar in one form or another, and 
the paper is written to help communication between experts in the various fields 
involved. The contents of Sections 2 and 7 and the numerical aspects of Section 5 
are, however, entirely new. 

2 Lower Bounds 

If Tt is a one-parameter semigroup with generator A, we define 
ujo = limsupfMogdlTtll), 

t— >+oo 

s = sup{Re(A) : A G Spec(A)}, 
Se = sup{Re(2;) : H-RzH > 
So = lim Se 

p = mm{uj : \\Tt\\ < e""* for all t > 0} 

where R{z) is the resolvent operator and e > 0. s and Sq are often called the 
spectral and pseudospectral abscissas respectively. An alternative characterization 
of p, sometimes called the logarithmic norm of A, is given in Lemma El One always 
has s < Sq < ujq < p, and each of these may be a strict inequality. In a Hilbert 
space uq = So, m Th. 5.1.11], so the value of ujq is deducible from knowledge of 
the pseudospectra (i.e. the resolvent norms). This identity is, however, not always 
valid in Banach spaces. 

The semigroup Tj (or its generator) is sometimes said to satisfy the weak stability 
principle if s = ujq, and the strong stability principle if there exists a constant M 
such that 

IITJ < Me'* 
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for all t > 0. Every diagonalizable matrix satisfies the strong stability principle, 
as does every operator in a Hilbert space which is similar to a normal operator. In 
Section 5 we will show that physically important self-adjoint operators need not 
satisfy the strong stability principle if they are considered with respect to a natural 
non-Hilbertian norm. 

In Example 1121 we show that ||Tj|| may oscillate rapidly with time. Because of this 
possibility we will not study the norm itself, but a regularization of it. Although 
our main application is to one-parameter semigroups, we work at a more general 
level to facilitate the discussions in the final section. We assume that B, T> are two 
Banach spaces and that Tt : V B is a. strongly continuous family of operators 
defined for t > 0, satisfying ||To|| = 1 and ||T(|| < Me*^* for some M,uj and all 
t > 0. We define N(t) to be the upper log-concave envelope of ||T(||. In other 
words uit) = log{N(t)) is defined to be the smallest concave function satisfying 
z/(t) > log(||Ti||) for all t > 0. It is immediate that N(t) is continuous for t > 0, 
and that 

1 = A^(0) < hm Ar(t). 

In many cases one may have N{t) = \\Tt\\, but we do not study this question, asking 
only for lower bounds on N{t) which are based on pseudospectral information. 

If k E H and we replace Tj by Tjc'^* then ||Ti|| is replaced by ||Ti||e'^* and N(t) is 
replaced by N{t)e'^^. We put k = —ujq or, equivalently, normalize our problem by 
assuming that ujq = 0. In the semigroup context this implies that Spec(A) C [z : 
Re{z) < 0}. It also implies that ||TJ > 1 for alH > by Th. 1.22]. If we 
define : V ^ B by 

POO 

RJ = / (r,/)e-^* dt 
Jo 

then II-R2II is uniformly bounded on {z : Re (z) > 7} for any 7 > 0, and the norm 
converges to as Re {z) — >• +00. In the semigroup context Rz is the resolvent of 
the generator A of the semigroup. 

The following lemma compares N{t) with the alternative regularization 

L{t) = sup{||T,|| : < s < t} 

of II Till, which was introduced by Trefethen, [21], and implemented in the package 
Eigtool by Wright, jSOl page 82], jSH. 

Lemma 1 If ujq = then 

\\Tt\\ < Lit) < Nit), 
for all t > 0. IfTt is a one-parameter semigroup then we also have 

N{t) < 

for all positive integers n and t > 0. 
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Proof The log-concavity of N{t) and the assumption that uq = imply that N{t) 
is a non-decreasing function of t. We conclude that ||TJ < L(t) < N(t). If Tj is a 
one-parameter semigroup we note that s — > L(t/?T,)^^""^/* is a log-concave function 
which dominates ||Ts|| for all s > 0, and which therefore also dominates N{s). 

In the following well-known lemma we put 

N'{0+) = \im_e-^{N{e) - N{0)} G [0, +oo]. 

Lemma 2 The constant p satisfies 

p = N\0+) > limsupt"^ {||Tt|| - 1} . 

IfTt is a one-parameter semigroup and B is a Hilbert space then 

p = sup{Re (z) : z e Num(A)} (1) 
where Num(A) is the numerical range of A. 

Proof If A^'(0+) < u then, since N(t) is log-concave, 

\\Tt\\ < N{t) < e'^* 

for all t > 0. The converse is similar. The second statement follows from the fact 
that, assuming A to be the generator of a one-parameter semigroup, A — uol is the 
generator of a contraction semigroup if and only if Num(A — uol) is contained in 
{z : Re (2) < 0}. 

We study the function N{t) via a transform, defined for all > by 

M{uj) = sup{||Tf||e-'^* : t > 0}. 

We see that up to a sign p{lj) = log(M(co')) is the Legendre transform of uit) 
(also called the conjugate function), and must be convex. It is also clear that 
M{uj) is a monotonic decreasing function of u which converges as u —* +00 to 
limsup^^g ||Ti||. Hence M{uj) > 1 for all u > 0. We also have 

N{t) = inf{M(cj)e'^* : < u < 00} (2) 

for all t > by the theory of the Legendre transform, or simple convexity argu- 
ments, fIE\ . 

In the semigroup context the constant c introduced below measures the deviation 
of the operator A from any generator of a contraction semigroup. 

Lemma 3 If a > 0, b eH and a\\Ra+ib\\ = c > 1 then 

M{u) > MM := I 5^ - < u; < r = a(l - 1/c) 

II otherwise. 
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Proof The formula 

roo 

Ra+^b = / T,e-("+*^)* dt (3) 
Jo 

implies that 

roo roc 

c/a < / N{t)e-''' dt < / M(cj)e'^*-''* dt 
Jo Jo 

for all u such that < uj < a. The estimate follows easily. 

This lemma is most useful when c is much larger than 1. If c = 1 then r = and 
the lemma reduces to M(a;) > 1 for all u; > 0. 

Theorem 4 // a||i?a+ib|| = c > 1 and r = a(l — 1/c) then 

N{t) > mm{e''\ c} 

for allt > 0. 
Proof This uses 

N{t) > inf{M(cu)e"* : > 0}, 

which follows from 

Trefethen and Wright give a related lower bound in (23 EO] , namely 

m> 



1 + (e»* - l)/c' 



Although these are lower bounds for different quantities and under slightly different 
conditions, the bound of Theorem H] is better in the following sense. The two sides 
of (jH) are asymptotically equal as t — 0+ and t ^ oo, but for intermediate t we 
have: 



Lemma 5 Let a > 0, t > and c > 1 then 

pit 

<min{e"(i-i/^)*,c}. (4) 



1 + (e«* - l)/c 

Proof Put s = at. There are two inequalities to prove for all s > 0. 

< c, 



6^ 



1 + (e^ - l)/c 

< o(i-i/c)s 
l + (e^-l)/c - 

After some algebraic manipulations, both are seen to be elementary. 

The above theorem provides a lower bound on N{t) from a single value of the 
resolvent norm. The well-known constants c(a), defined for a > by 

c(a) = a sup{ II i?a+if> II : b G R}, 



5 



are immediately calculable from the pseudospectra. it follows from Q and uq = 
that c{a) remains bounded as a — +00. The transform c(-) defined below is easily 
calculated from c(-). 

Corollary 6 Under the above assumptions one has 

N{t)>d{t):= sup (min{e"(")*, c(a)}}. (5) 

{a:c(a)>l} ^ ^ 

where 

r(a) = a(l — l/c(a)). 

Theorem 7 // Tf is a one-parameter semigroup and sq = ujq = then c(a) > 1 
for all a > 0. 



Proof If c(a) < 1 then by using the resolvent expansion one obtains 

II 7-1 II c(oj] , I I / \ / \— 1 

-Ka+i6+2 < (1 - \z\c{a) a) 

a 

for all \z\ < a/c{a). This implies that sq < 0. 

Examples show that c(a) is often a decreasing function of a, but this is not true 
in Example El below. The supremum in (jSJ need only be taken over those a at 
which c(a) is decreasing. The following transform of c(-) may sometimes be easier 
to compute than c(-). 

Lemma 8 //c(-) > 1 is a monotonia decreasing function then N{t) > c{t) for all 
t > 0, where c(-) is the function inverse to 

tic) ■ 



a(c)(l - 1/c) 



and c — > a(c) is the function inverse to a ^ c(a). The function c{t) is defined for 
alio <t < 00. 



Proof We only consider the case in which c(-) is differentiable with a negative 
derivative at each point. The graph of t — ^ N{t) li^s above each of the points 
(t(c),c) by Theorem m Putting g{c) = log(c)/(l — 1/c) a direct calculation shows 
that g'{c) > for all c > 1. Its definition implies that a'(c) < for all c. Hence 

,// N aic)g'{c) - a'{c)g{c) 
'^'^ = >° 

The domain of c(-) is the same as the range of t(-), and one may show that t{c{a)) — »• 
as a — > 00, while t(c(a)) ^ cxd as a ^ 0. 

Although much recent progress has been made, the numerical computation of the 
pseudospectra is still relatively expensive. All of the examples of one-parameter 
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semigroups in the next section are positivity-preserving, in the sense that / > 
imphes Tj/ > for all t > 0. In this situation the evaluation of c(a) is particularly 
simple. The following is only one of many special properties of positivity-preserving 
semigroups to be found in \t^ . 

Lemma 9 Let Tt be a positivity-preserving one-parameter semigroup acting in 
LP{X, dx) for some 1 < p < oo. If = then 

\\Ra+ib\\ ^ ||-Ra|| 

for all a > and 6 G R. Hence c{a) = a\\Ra\\. 

Proof Let f e and g e L" = {W)* , where l/p + 1/g = 1. Then 

\{RaM^f.g)\ = I / {TJ,g)e-^^^"'^Ut\ 
Jo 

/•CO 

< / \{TJ,g)\e-''' dt 
Jo 

/•oo 

< mi\g\)e~'^Ut 

= {Ra\f\,\9\) 

< ll^allll/llplkll,. 

By letting / and g vary we obtain the statement of the lemma. (The inequality 
\{Xf,g)\ < \g\) for all positivity-preserving operators X may be proved by 

considering first the case in which f,g take only a finite number of values.) 



3 A Direct Method 

The direct calculation of ||Tt|| for t > is not straightforward for very large ma- 
trices, i.e in dimensions of order 10^, particularly when using the norm: even if 
A is sparse, e^^ is usually a full matrix. If the generator A of the semigroup has 
enough eigenvalues the following method may be useful. Let {/r}"=i be a linearly 
independent set of vectors in B, and suppose that Afr = Xrfr for 1 < r < n. Let 
C denote the linear span of {/i, /„} and let Tc,t denote the restriction of Tt to 
C It is clear that 

II^J > \\Tc,t\\ 

for all if: > 0. If £ is large enough one might hope that this is a reasonably good 
lower bound. If A has a large number of eigenvalues, then one might choose some 
of them to carry out the above computation after inspecting the pseudospectra of 
A. 

The operator Tjr^t must be distinguished from PnTfPn, where P„ is the spectral 
projection of A associated with the set of eigenvalues {Ai,...,A„}. If n = 1 the 
norm of the first operator is |e~^i*| while the norm of the second is ||Pi|| |e~^i*|. 
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We will see in Table 4 that the norm of Pi may be very large. Unfortunately the 
norm of PnTtP^ is much easier to compute than that of T^^t in the context, using 
Matlab's eigs and norm(-,l) routines, and it is easy to confuse the two. 

The following standard result is included for completeness. 

Lemma 10 Under the above assumptions, suppose also that the linear span of 
{/r}^i is dense in B. Let T^^t denote the restriction ofTf to Cn — ■■■,/n}- 
Then 

lim ||T„,t|| = ||Tt|| 

for all t > 0. If t ^ ||Tt|| is continuous on [a, 6] then the limit is locally uniform 
with respect to t on that interval. 



Proof Given e > and t > there exists f E B such that ||/|| = 1 and ||Tt/|| > 
\\Tt\\ — e. By the assumed density property, we may assume that f E Cn for some 
n. This immediately yields ||Tt|| > ||T„_t|| > \\Tt\\ — e. 

The final statement is a general property of any pointwise, monotonically conver- 
gent sequence of continuous functions to a continuous limit. 

Clearly this lemma is of limited use in the absence of any information about the rate 
of convergence. If B is a Hilbert space, the norms of the approximating semigroups 
may be evaluated by the following standard result. We know of no analogue of 
this lemma for subspaces of Banach spaces. The problem is that the unit balls of 
subspaces of may have very complicated shapes, which makes operator norms 
difficult to compute. For example the unit ball of a generic, real, two-dimensional 
subspace of n} is a polygon with 2n sides, and higher dimensional subspaces 

are even more complicated. 

Lemma 11 If B^ g = {fs, fr) (^nd D^- g^t = ©^''^^r.s for 1 < r, s < n, then 

where the norm on the RHS is the operator norm, C" being provided with its 
standard inner product. 



Proof The n x n matrix B is readily seen to be self-adjoint and positive. If S : 
C" Cn is defined by 

n 

Sa = ^Prfr 
r=l 

where (5 = B~^l'^a., then S is unitary and 
This yields the statement of the lemma. 
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4 Exactly Soluble Examples 



Example 12 Let Tt be the positivity-preserving, one-parameter semigroup acting 
on L^(R+) with generator 

Af{x) = f\x) + v{x)f{x) 
where v is any real- valued, bounded measurable function on R"^. Explicitly 

T,f{x) = ''-^^f{x + t) (6) 

for all f E and alH > 0, where 

a{x) = exp |y v{s) ds 

The function a is continuous and satisfies 

e"ll''ll°°*a(x) < a{x + t) < ell''l'°°*a(x) 
for all x,t. hence \\Tt\\ < ell''ll°=* for all t >0. 

The precise behaviour of ||Tt|| depends on the choice of v, or of a, and there is a 
wide variety of possibilities. For example if c > 1 and b > then the choice 

a{x) = l + {c-l)sm^{7cx/2b) (7) 

leads to ||T2„b|| = 1 and ||T'(2„+i)fe|| = c for all positive integers n. In the case ((7j), 
the regularizations N{t) and L{t) are not equal, but both are equal to c for t > b. 

If c > and < 7 < 1 then the unbounded potential v{x) = c(l — 'y)x^'^ corre- 
sponds to the choice 

a{x) = expjcx"^"'*'}. 

Instead of deciding the precise domain of the generator A, we define the one- 
parameter semigroup Tf directly by ©, and observe that 

N{t) = \\Tt\\ = exp{ct^"^} 

for all t > 0. If c is large and 7 is close to 1, the semigroup norm grows rapidly 
for small t, before becoming almost stationary. The behaviour of \\Ttf\\ as t 00 
depends upon the choice of /, but for any / with compact support Ttf = for all 
large enough t. On the other hand ||Tt/|| cannot be a bounded function of t for all 
/ G L^(R"'"), because of the uniform boundedness theorem. 

For this unbounded potential v, every z with Re {z) < is an eigenvalue, the 
corresponding eigenvector being 

f{x) = exp |zx — c(l — '-f)x^~^^ . 
9 



Hence 

Spec(A) = {z : Re{z) < 0}. 

On the other hand p = +00, and Num(A), which is always a convex set, must 
equal the entire complex plane by Lemma |2l 

Example 13 If we put 

^"[00 

acting in with the Euclidean norm, and |A| = r > then 

1 ^1 + l/4r2 



so 



2r2 



1 



c(a) = — + ^1 + l/4a2 



for all a > 0. We also have 



\Tt\\=t/2 + Jl + tyA 



which is log-concave, so \\Tt\\ = N{t) for all t > 0. The choice a = 2 provides a 
fairly good lower bound on N{t) for < t < 0.5. As a gets smaller we get a better 
lower bound on N(t) for large t > 0, while as a gets bigger we get a better lower 
bound for small t > 0. 

Example 14 Let A be the n x n Jordan matrix 



1 if j = i + 1 
otherwise 



acting in C" with the norm. Then 

Tt = I + At + + ... + A"-H"-V(n - 1)! 

and 

llTtll = l + t + e/2\ + ... + f'-^/in - 1)! 
for all t > 0. A direct calculation shows that 

\\Tt\\=N{t) = L{t) 

for all t > 0, and that p = 1. 

Direct calculations are not so easy for this example if one uses the norm. 
However, in this case it follows from (Q) that p equals the largest eigenvalue of 
B = {A + A*)/2. Since the set of eigenvalues is {cos(r7r/(n + 1))}"^;^, it follows 
that 

p = cos(7r/ {n + 1)) < 1. 
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Example 15 Given 7 > 0, let 



-7 


1 








-7 















act in C'^ with the Euclidean norm. We have 

\\Tt\\ =max|l,e-^*|t/2 + yTTtV4||. 

If < 7 < 1 then this is not log-concave, and for 7 close to 0, it increases linearly 
in t for a long time, before eventually dropping to 1. The functions L{t) and N{t) 
are equal, and they are constant for large enough t > 0. In this example 

c(a) = max < 1, — H — Jl H — ; — ■ I . 

^ ' \ '2(a + 7)2 a + 7V 4(a + 7)2/ 

This equals 1 for small a > 0, and converges to 1 as a ^ 00, but it is not a 
monotonic decreasing function of a. 



5 Schrodinger Semigroups 

Semigroups with generators of the form A = —H = /\ — V have been extensively 
studied, and provide a fascinating insight into the importance of the Banach space 
on which they are chosen to act. 

If one assumes that the potential (multiplication operator ) V : — R lies in 
the so-called Kato class, then the self-adjoint operator H = —A + V may be inter- 
preted as a quadratic form sum in ^^(R^), and the one-parameter 'Schrodinger 
semigroup' {e~'^*}t>o on may be extended consistently to all of the spaces, 
1 <p < 00, 1^. 

If H is interpreted as a quantum-mechanical Hamiltonian, then there are good 
reasons for being interested primarily in the choice p = 2. We show in the next 
section that the time- dependent Schrodinger equation f'{t) = —iHf(t) is only 
soluble in for p = 2, but in addition the use of the norm is fundamental to 
the probabilistic interpretation of quantum mechanics. In this context Schrodinger 
semigroups are only of technical interest; they enable one to investigate a variety 
of spectral questions very efficiently. 

When studied in the spectral theorem yields the strong stability condition 

||e-^*|| =e-"* 

where 

A = min{Spec(if)}. 
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If the potential V depends upon a parameter c, then one often has A(c) = for 
some range of values of c, with transitions to A(c) < at certain critical values of 
c. These critical values describe the sharp emergence of instability. 

Schrodinger semigroups also have direct physical significance in problems involving 
diffusion, but in this context the equation should be studied in L^(R^). The point 
here is that the semigroup e~^* is positivity-preserving, and ft = e~^^ describes 
the distribution of some continuous quantity in at time t > given its initial 
distribution /. Assuming / > 0, the total amount of the quantity at time t is given 
by 

/ /(t,x)d^x=||/i||i. 

It is known that, in the technical context described above, the spectrum of Hp (the 
operator H considered as acting in W) is independent of p, The operators 

^-Ht g^j^g known to have positive 'heat' kernels K{t,x,y), [SH!, and 



|e-^i*|| = sup / K{t,x,y) d^x. 



We will see that these integrals of the heat kernel are not determined by the 
spectral properties of Hi. We start by showing that the value of the constant 
p may be entirely different in the and contexts. The conditions on the 
potential V in the following theorem can clearly be weakened, and we refer to P] 
for a comprehensive treatment of the problem. 

Theorem 16 Let Hi = —A + V , where V is continuous and bounded below, with 

c = mf{V{x) : X E R^}. 

Then c = —p. 

Proof The inequality p < — c, or equivalently 

||e--f^i*|| < e""* for alH > 
may be proved by the use of functional integration or the Trotter product formula. 



Conversely, let e > and let |a; — a| < 6 imply c < V{x) < c + e. Let / G C^{{x 
\x — a\ < 6}) be non-negative with ||/||i = 1. Then 

= limt-Mm/,1) -(/,!)} 

= -{Hif,l) 
= {Af-Vf,l) 

= -{vf,i) 

> -(c + 5)(/,l) 

= —c — e. 
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This implies that p > —c. 

Corollary 17 // Hi = —A + V where V is not bounded below, then p = oo, 
whatever the spectral properties of Hi . 

The above results show that the short time semigroup growth properties do 
not depend only upon whether the spectrum is non-negative. We cannot give 
a complete analysis of the long time behaviour, since the requisite theorems do 
not exist, but discuss a typical case below. Our main purpose is to emphasize 
that one may have a failure of the strong stability principle for such semigroups. 
Generalizations of this example have been studied in considerable detail by Murata, 
[TKt I16j. and by Davies and Simon, jBI. The most general results which we know 
about are by Zhang, 



< c< ^ ^. 

4 



Example 18 Let > 3 and let 



so that 




2 + 

Now consider the operator Hp = —A + V acting in L^(R^), where the bounded, 
strongly subcritical potential V is defined by 



V{x) 



-c\x\ 2 if |x| > 1 
otherwise. 



It is known that the operator Hp has spectrum [0, oo) for all 1 < p < oo, and that 
c = (A^ — 2)2/4 is a critical value for the emergence of a negative eigenvalue, [T^ I5] 



The operator Hp possesses a zero energy resonance rj given by 
< r]{x) = \ 



if |x| > 1 
1 — P otherwise. 



where 

</3 = — < 1. 

The operator —Hp generates a positivity-preserving one-parameter semigroup act- 
ing in LP(R^) for all 1 < p < oo, and for p = 2 it is a self-adjoint contraction 
semigroup. On the other hand it is proved in jH| Th. 14] that for any ai, (72 
satisfying < (Ti < q;^/2 < < cxd there exist positive constants Ci, C2 such that 

ciil + ty^ < ||e-^^*|| < C2(l + t)"^ (9) 

for all t > 0, the norm being the operator norm in Li(R^). We conclude that 
s = Sq = ujq = for this example, whether the operator is considered to act in 
LHR"") or L^fR^). 
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The above example exhibits polynomial growth of the operator norm as i — > oo. 
It exhibits the weak, but not the strong, stability property. 

Theorem 19 For every 7 > there exists a Schrddinger semigroup e~^p^ acting 
in i/(R-^) for all 1 < p < 00 such that 

ci(l + -fHy < ||e-^i*|| < C2(l + -fHy^ 

for all t > 0, even though K2 ^K*>0 in L^{R^). 

Proof The operator is given by Kp = —A + V^, where 




—c\x\ ^ if > I/7 
otherwise. 



The bounds are proved by reducing to the case 7 = 1 by using the scaling trans- 
formation {U^f){x) = 7^/V(7^)- 

By exploiting the rotational invariance, it is easily seen that the above example 
is associated with a similar example on the half-line. However, the transference 
procedure is different for the L} and norms. 

Lemma 20 Let the potential V he rotationally invariant and hounded helow on 
R^. Then the self-adjoint operator H — —A -|- V , defined as a quadratic form 

sum, is hounded helow, and the one-parameter semigroup Tt defined for i > 
hy Tt = e^^* acts consistently on L*'(R^) for all 1 < p < 00 and commutes 
with rotations. If we identify the rotationally invariant suhspace o/L^(R^) with 
L^((0, 00), dr) in the usual way, then the restriction of H2 to this suhspace is given 
hy 

(At iL'T^ 

subject to Dirichlet boundary conditions atr = 0. On the other hand if we identify 
the rotationally invariant suhspace of L^{Il^) with -^^^((O, 00), dr) in the usual way, 
then the restriction of Hi to this suhspace is given hy 

^i/W = + - Wir)/ry + V{r)f{r) 

subject to Dirichlet boundary conditions at r — 0. 

Proof The operator H acts on the space -£'^((0, 00), r^~^ dr) of rotationally invari- 
ant functions according to the formula 

^2/(r) = p-^f } + F(r)/(r). 

We transfer H2 to L^{{0, 00), dr) by means of the unitary map Uf{r) — r^^~^^^'^f(r), 
obtaining the stated formula for L2 = UH2U~^. 
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The operator Hi acts on the space L^((0, oo), r^~^ dr) of rotationally invariant 
functions according to the same formula as for H2- We transfer Hi to ((0, oo) , dr) 
by means of the isometric map V f{r) = r^~^/(r), obtaining the stated formula 
for Li = VHiV-\ 

There are two ways of seeing that one should impose Dirichlet boundary conditions 
at r = 0. If one calculates the heat kernels using functional integration, the relevant 
fact is that the probability of Brownian motion passing through the origin in 
for > 2 is zero, j2Sl- Alternatively, subject to minimal regularity conditions on 
/ at the origin, we see from their definitions that Uf{0) = Vf{0) = 0, at least if 
> 3. 

There are several ways of discretizing the operator Li. One obtains a discretiza- 
tion which has real eigenvalues and generates a positivity-preserving semigroup by 
starting from the formula 

LJ{r) = r(^-i)/2L2{r-(^-i)/V(r)}. 

The last part of the following lemma will be used when carrying out numerical 
calculations below. 

Lemma 21 Let M2 be a self-adjoint n x n matrix with non-positive off-diagonal 
entries, and let D he a diagonal nxn matrix with positive entries. Then the matrix 

Ml = DM2D-^ 

has the same, real, spectrum as M2, ande~^^^* is positivity-preserving for allt > 0. 
// also M^l > then e~^'^^'^ is a contraction semigroup on C" provided with the 
norm. If A is an eigenvalue of M2 with multiplicity 1 and f ^ is a corresponding 
eigenvector, then the spectral projection Px of Mi corresponding to the eigenvalue 
A has norm 

||p||_ P/l|l||p-^)7l|oo 

" (/,/) 

calculated with respect to the norm o/C". 



(10) 



Proof See jSl Th. 7.14] or the proof of Lemma for the positivity- preservation. 
The second statement is also classical, but we include a proof for completeness. 
Since the coefficients of the matrix e"*^^* are non-negative, we have 



■n 

|e-*'^Vl|i-||/l|i = E{l(e"''^V).|-IM} 

E{(e-^^*I/I).-|M} 



n 

< 

r=l 

(e-"'^*l/|- 1/1,1) 

\Mie-^'^V\,l) ds 











\e-''-V\,M*l) ds 



< 0. 
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The expression for ||Pa|| is obtained from the formula 



if J) 



Df. 



Example 22 We describe a discretization of the operator Li, with the critical 
value of the parameter c in (jHJ, namely c = [N — 2)^/4, and with = 3, acting in 
the space C" of finite sequences. We put 

r (2-i;i)/i-/2 ifr = l 

(M2/), = (2 - Vr)fr - fr^l -fr+1 if 2 < r < 72 - 1 
[ (2 - Vn)fn - fn-1 Hv = n 

We choose 

^'^ ~ j 2 - s;\ ~Sr if 2 < r < n 
where = (1 + l/r)"*^/^. We note that 

as r — > CX3. We finally put Mi = DM2D^^ where Dr^s = r6r,s for all r, s. 

Theorem 23 The matrix M2 is non-negative and self-adjoint. The operators 
Q-Mit Qn positivity-preserving for all t > 0, and their eigenvalues Xr^t M 
satisfy < \.^t < 1- If we replace Vr by in the above definitions, then e~^i* is a 
one-parameter contraction semigroup on C"" provided with the norm. 

Proof The self-adjointness of M2 is evident. The fact that M2 is non-negative 
depends upon a discrete analogue of the Hardy inequality. There is a substantial 
literature on discrete analogues of differential inequalities, but we can prove the 
result which we need very quickly. The relevant quadratic form is 



n 

Q{a) = |aip + |a„+ip + ^ ||ar - flr-il 



Uip I Q/f 



|2 I 



r=l 
n 



y^(2 — t>r)|ar|^ — ^{a^Cr-l + Ur-lUr} 
r=l r=2 



El 1/2 -1/2 |2 , I 1/2 

\Sj._iar-l — Sr_i Qrl + |S„ Qn} 



r=2 
> 0. 



Since Ml and M2 are similar, the comments about the eigenvalues of e~^^* follow 
immediately. The fact that e~^i* is positivity preserving for t > follows using 
Lemma I?T1 as does the final statement of the theorem. 

In spite of the above. Theorem ITIIl suggests that the norm of e~*^^*, considered as 
an operator on C" provided with the norm, should grow with t. Table 1 shows 
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the results of testing this numerically using Matlab. Our computations used the 
formula 

and exploited the self-adjointness of M2 when calculating the exponential. One 
can use this formula directly to obtain the bound 

||e-Mit|| < nV2|p|| ||e-M,t|| < ^3/2 

for all t > 0, using the fact that ||/||2 < ||/||i < nV2||/||2 for all / G C" and 
M2 — M| > 0. However, this provides no insight into the hmiting behaviour as 
n 00. 

Table 1. Values of ||e^^^^*|| for various n 
t n = 100 n = 200 n = 300 






1 


1 


1 


100 


4.059 


4.059 


4.059 


200 


4.824 


4.824 


4.824 


300 


5.333 


5.337 


5.337 


400 


5.701 


5.735 


5.735 


500 


5.945 


6.063 


6.063 


600 


6.071 


6.346 


6.346 


700 


6.095 


6.595 


6.595 


800 


6.036 


6.818 


6.818 


900 


5.914 


7.022 


7.022 


1000 


5.747 


7.208 


7.209 



For n = 300 the maximum value of the norm occurs for t ~ 6000. While the 
increase may not appear very rapid, it should be noted that we have assumed 
a unit separation of the points on Z"^, so the implied time scale is very long by 
comparison with that of the corresponding differential operator. Table 2 shows 
how the maximum value of the norm as t varies depends upon the value of n. 

Table 2. maxj>o He^*"^^*]! as a function of n. 

n norm max 

^0 433 

100 6.10 

150 7.46 

200 8.60 

250 9.61 

300 10.53 

Since every eigenvalue of Mi is positive we must have 

lim ||e-^^*|| = 

t— >oo 
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but, still using the norm, if n = 300 the inequality ||e ^^*|| < 1 only holds for 
t > 4.6 X 10^ 

Table 1 suggests the existence of a limit as n ^ oo, and this is proved below. We 
identify with the subspace of /^(Z"*") consisting of sequences with support in 
{1, ...,n}. We also identify any n x n matrix X with the operator X on /^(Z+) 
defined by 

(Xf) =i ^-1^-'^^^ ifl<r<n 
^ [0 otherwise. 

We finally exhibit the n-dependence of the various operators explicitly. 

Lemma 24 There exists a bounded operator Mi ^^o on 1^(7a^) to which Mi^„ con- 
verge strongly as n ^ oo. For every t > the operators e~^i'"* increase monoton- 
ically to e~^^i'°°*, and 

lim ||e-^'^^'"*|| = ||e-*^i-°°*||. (11) 



Proof The limit operator is given by 

(Mi,oo/)r = 



{2-v{)h-\h ifr = l 

(2 - i;,)/. - ^/.-i - ifr>2 



l,oo 



and is evidently bounded on / (Z"*^). The strong convergence of Mi^„ to M; 
implies the strong convergence of the semigroup operators. We also have 

oo 

for all / G /^(Z+), where K„,{t,r,s) > is the transition 'probability' for a jump 
process which is killed if it moves outside {1, n} and grows at the rate Vr at each 
r such that 1 < r < n. It follows on probabilistic grounds that n Kn{t,r,s) is 
monotonic increasing with 

\un^Kn(t,r, s) = K^{t,r,s). 

This implies (jlll) . 

The formula © with N = 3 suggests that for our example one should have 

||e-Mi,^t|| _ ^^1/4 

as t — i> oo. For finite n this can only happen for t in the transitory growth interval. 
Numerical calculations confirm this. If n = 200 one has 

(2.69t)i/' < ||e-^^i*|| < (2.71t)^/' 

for all t satisfying 200 < t < 1200. If n = 300 the same holds for 200 < t < 2500. 

We may also investigate the resolvent norms in the context, or more specifically 
the function c(a) = a||-Ra||; see Lemma El The eigenvalues of Mi are all positive. 
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so ojq 7^ 0, and we must have lima^o+ c(a) = 0. However, the smallest eigenvalue 
converges to as n ^ oo, so c(a) may be quite large even for small positive a. 
The data in Table 3 were obtained for the case n = 300, putting a = 2"™" and 
stopping at the value of m for which c(a) takes its maximum value. For n = 300 
the smallest eigenvalue of Mi is 6.38 x 10~^ and the largest is 4.00. 

Table 3. Dependence of c(a) on a = 2~™ for n = 300. 

m c(a) 
~ 1.50 

2 1.72 

4 2.36 

6 3.30 

8 4.65 
10 6.56 
12 8.45 

For n = 1000 the smallest eigenvalue of Mi is 5.772 x 10^^ and the largest is 4.00. 
The largest value of c(a) for a of the above form occurs for a = 2^^^ and is 15.50. 

We finally tabulate how the smallest eigenvalue A of Mi depends upon n, with the 
values of the norm of the corresponding spectral projection Pa, computed using 
(fTUI) . The fact that ||Pa|| grows like n^^^ as n increases was expected on the basis 
of replacing / in fjlOj) by the exact zero energy resonance gr = r^^"^ of the operator 
Ml acting in L^(Z+). 

Table 4. Dependence of A and ||Pa|| on n. 



n 


A 


II^aII 


\\Px\\/n'/' 


100 


5.669 X 10"^ 


11.178 


1.1178 


200 


1.4314 X 10"^ 


15.772 


1.1152 


400 


3.597 X 10-^ 


22.278 


1.1139 


600 


1.601 X 10"^ 


27.274 


1.1134 


800 


9.014 X 10-*^ 


31.487 


1.1132 


1000 


5.772 X lO-'^ 


35.199 


1.1131 



Example 25 In the above study we focussed on the case = 3, but the difference 
between the and theories becomes even more dramatic for larger values of A^. 
The only change needed in our discrete example with the critical value of c in (jH)), 
namely c = (A^ — 2)^/4, is to redefine D by D^^s = r^^~^'^^'^5r,s for all r, s. For = 6 
the bounds Q then suggest that ||e~^^i*|| ~ t as t ^ oo. Numerical calculations 
yield 

4.00 1 < ||e"^^i*|| < 4.02 1 
for all t satisfying 100 < t < 2000, when n = 300. 
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Conjecture Let > 2, let {Df)r = r^^-^^^fr for all r > 1, and let 

{M2,oof)r = 



{2-vr)h-h ifr = l 

(2 - Vr)fr - fr-1 " /r+1 if T > 2. 



Then Mi^oo = DM2.00D ^ is a bounded operator on /^(Z"*") with non-negative real 
spectrum, and there exists a positive constant c such that 

hm t-(^-2)/4||e-*^^-*|| =c. 

t— >oo 



6 Absence of Upper Bounds 

In finite dimensions it is also possible to obtain upper bounds on semigroup norms 
from spectral or pseudospectral information, but the results deteriorate as the 
dimension increases, El IH 123 • It is therefore not surprising that no such bounds 
can be obtained in a general Banach space setting. In this section we describe 
physically important examples to show that this difficulty cannot be evaded. 

The converse part of the following theorem is a classical result of Hille and Yosida, 
and has frequently been used to pass from resolvent bounds or from the dissipative 
property to a one-parameter semigroup, jSl Cor. 2.22]. The smallest possible 
constant c in (fT^ is often called the Kreiss constant by numerical analysts, by 
analogy with the constant of the Kreiss matrix theorem. 

Theorem 26 IfTt is a one-parameter semigroup satisfying \\Tt\\ < c for allt>0 
then its generator A satisfies 

Spec(A) C {A : Re (A) < 0} 

and 

for all A such that Re (A) > 0. The converse implication holds if c = 1. 

There are many important examples in which one does not have c = 1. The 
following is typical of semigroups whose generator is an elliptic operator of order 
greater than 2, and is treated in detail in [Hj. 

Example 27 Let Tt act in L^(R") for t > according to the formula 

Ttf{x) = h * fix) 
where * denotes convolution and 

4(0 = e-'«'^* 
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Formally speaking Tt = e^^ where A = — A^. It is immediate that kt lies in 
Schwartz space for every t > 0, and hence that convolution by kt defines a bounded 
operator on L^. kt is not a positive function on R", and if we put c„ = then 
> 1 is independent of t by scaling and 

\\Tt\\ = Cn 

for all t > 0. For n = 1 we have ci ~ 1.2367. 



The following more general theorem implies that p = +oo for all one-parameter 
semigroup whose generator is elliptic of order greater than 2, jl4|. 

Theorem 28 Let Q be a region in and let A be an elliptic operator of order 
greater than two whose domain contains C^{Q). If A generates a one-parameter 
semigroup Tt on U'{VL) and p ^ 2 then Tt cannot be a contraction semigroup. 

In spite of its great value, we emphasize that the Hille-Yosida theorem is numer- 
ically fragile. An estimate which differs from that required by an unmeasurably 
small amount does not imply the existence of a corresponding one-parameter semi- 
group. We conjecture that an example with similar properties can be constructed 
in Hilbert space. 

Theorem 29 For every e > there exists a reflexive Banach space B and a closed 
densely defined operator A on B such that 

(i) Spec(A) C iR, 

(ii) ||(AJ- < (1 +e)/|Re(A)| forallXiiH, 
(Hi) A is not the generator of a one-parameter semigroup. 

Proof Given 1 < p < 2, we define the operator A on L^(R) by 



Af{x) 



dx^ 



As initial domain we choose Schwartz space S, which is dense in L^(R). The 
closure of A, which we denote by the same symbol, has resolvent operators given 
hy R\f = gx * f , where * denotes convolution and 

for all A ^ m. If p = 2 the unitarity of the Fourier transform implies that 
||-Ra|| < |Re (A)!^"*^. For p = 1, however, 

II-RaII = I^aIIli- 

Assuming for definiteness that Re (A) > the explicit formula for gx yields 

1 



i f°° r 

|i?A|| = exp [-|x|Re{(2A)i/2} 



dx. 
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Putting A = re*^ where r > and — 7r/2 < 9 < 7r/2, we get 



rcos(0/2 + 7r/4) " |Re(A)r 

Interpolation then imphes that if 1 < p < 2 and l/p = 7 + (l — 7)/2 then 

27 

"^^1' - WW 

By taking p close enough to 2 we achieve the condition {ii). 

The operators Tt are given for t 7^ by Ttf = kt * f where * denotes convolution 
and 

kt{x) = {Anity^^'^ expi-x"^ /Ait}. 

It follows from the formula for the operator norm on i^^(R) that Tt are not bounded 
operators on L^(R) for any t ^ 0. Suppose next that 1 < p < 2 and that a 
semigroup Tt on L^(R) with generator A does exist; we will derive a contradiction 
by an argument which goes back at least forty years. If / G 5 and ft & S is defined 
for alH G R by 

then ft is differentiable with respect to the Schwartz space topology, and therefore 
with respect to the norm topology, with derivative Aft. It follows by [Si Th. 1.7] 
that ft = Ttf. Now assume that a > and /(O = so that ft{^) = e-('*+**)«'. 

Explicit calculations of ft and / yield 

p 



Hence 



(47ra)i/2p-i/2p-i/2p 

||/,||^ = (47r)l/2p-V2p-l/2p^-l/2p(^2^^2)l/2p-l/4_ 
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\Tt\\>fjf^ = {l + tya'f-^y'^. 
Il/llp 

But this diverges to 00 as a 0, so cannot exist as a bounded operator for any 
t ^ 0. 

The above theorem implies that one cannot expect to derive upper bounds on 
semigroup norms from numerical resolvent norm estimates, i.e. from pseudospec- 
tra, in infinite-dimensional contexts. The Miyadera-Hille-Yosida-Phillips theorem 
provides a general connection between resolvent and semigroup bounds, jSl The- 
orem 2.21]. However, it involves obtaining bounds on all powers of the resolvent, 
and is rarely useful. 



7 Special Initial States 

It might be hoped that the pathologies described above are a result of applying 
the semigroup to 'untypical, badly behaved' initial vectors /, and that they would 
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disappear if / is restricted in an appropriate manner. In this section we investigate 
the consequences of assuming that / hes in the domain of A, so that t Tff 
satisfies the Cauchy problem in the classical sense. This amounts to studying the 
behaviour of Tj regarded as an operator from T> := Dom{A) to B, where the former 
space is provided with the natural Banach space norm 

Ill/Ill = MII/f + P/f)'/^ 

and k is chosen so that the embedding operator from V to B has norm 1. Almost 
equivalently one can study the regularized operators Tt = hTfRa on B for some 
a ^ Spec(y4), where h = The lower bounds of Section 2 are apphcable to 

either of these families of operators, the appropriate 'resolvent' operators in the 
second case being 

Rz = hRzRa = {Rz — Ra) ■ (13) 

a — z 

In the first case the resolvent operators are unchanged, but the values of their 
norms change. It follows immediately from (fT!?|l that h — > ||i?a+jf,|| is bounded if 
and only if 6 ^ (1 + |6|)||_Ra+jf,|| is bounded. 

Lemma 30 The operators Tt depend norm continuously on t and satisfy the short 
time growth condition 

\\ft\\<l+tL{t){a+\\RaV). (14) 

Proof The norm continuity oit-^Tt follows from the formula 

{ft-t)f = hjy,{ARa)fdx 

= h T^ {aRa - 1) / dx, 

proved using |^ Lemma 1.2]. This implies ()14|) by putting s = 0. The following 
theorem is similar to a result in J2] , and both are implied by Theorem |221 below. 

Theorem 31 Assuming that A is unbounded, one has 

Spec(Tt) = {0} U {he^'ia - X)-^ : A G Spec(A)} 

for all t > and a > cjq. 

Proof We normalize the problem by putting A = A — where < 7 < a, 
a = a — 7 and Tt = e~'^^Tt, so that 

TtR, = e^'ftka. 

The semigroup Tt is uniformly bounded since cjq = 1^0 ~ 7 < 0. Moreover 

TtRa = / f{s)T, ds 
Jo 
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where 

^, , f if < s < t 



e 



-ais-t) [{ s> t. 



Since / G L^{0, oo), the stated result is imphed by our next, more general, theorem, 
which appears to be new. 

Theorem 32 Let Tf be a uniformly bounded one-parameter semigroup acting on 
B, with an unbounded generator A. Let f G -^^^(0, oo) and 

POO 

Xf = / f{t)T, dt, 

<j 

where the integral converges strongly in C{B). Put 

f{z) = / /(t)e^* dt 
Jo 

for all z satisfying Re (z) < 0. Then 

Spec(Xj) = {0} U {/(A) : A G Spec(A)}. 



Proof We follow the approach of Th. 2.15]. Let ^ be a maximal abelian 
subalgebra of C{B) which contains Tt for alH > and Ra for all a ^ Spec(y4). Let 
M denote its maximal ideal space of A and " the Gelfand transform. Then A is 
closed under the taking of inverses and strong operator limits. Hence Xf G A and 

Spec(D) = {D{m) : m e M} 

for all D e A. 

If a, 6 ^ Spec(y4) then the identity 

Ra{m) - Rb{m) = (6 - a)Ra{m)Rb{m) (15) 
implies that the closed set 

N = {meM : Raim) = 0} 

is independent of the choice of a. Since A is unbounded must be non-empty. If 
m G M\N then 

Raim) G Spec(/?J\{0} = (a - X^)-' 

for some Am G Spec (A). A second application of (fTS|) implies that Am does not 
depend upon a. The definition of the topology of M implies that A : M\N — > 
Spec (A) is continuous. 

Let V denote the set of all functions / : (0, oo) — >^ C of the form 

n 

/(t) = ^a,e-^^* 

r=l 



24 



where Re > for all r. For such a function 

n 

Xf = Y.^rRl3r- 
r=l 

Therefore 

Spec(X/) = [Xf{m) : m G m} 

= {0} U |x^^a,i?^,,(m) : m G M\iv| 

= {0}u(x]«r(/5r-A™)-^:mGM\iv| 



>r=l 
' n 



= {0}u|^^^a,(/5, - A)~i : A G Spec(A)| 

= {0} U |^°° /(t)e^* dt : A G Spec(A)| 
= {0}U{/(A): AGSpec(A)}. 

Finally let / be a general element of L^(0,oo). There exists a sequence /„ G P 
which converges in norm to /, and this implies that Xf^ converges in norm to 
Xf, and that /„ converges uniformly to /. Hence 

Spec(X/) = JimSpec(X/J 

= {OjUjim {A(A):AGSpec(A)} 
= {0}u{/(A):AGSpec(A)}. 



In this final step we used the fact that {0} U {/(A) : A G Spec(74)} is a closed set. 
This is because Spec(74) is a closed subset of {2; G C : Re {z) < 0}, and f{z) 
as 1^1 — 00 within this set. 

In spite of Theorem El Wrobel's modification of the example of Zabczyk shows 
that the long time growth properties of Ti cannot be deduced from its spectral 
behaviour, [3^1 Ex. 4.1]; a special case is described in [2, Ex. 5.1.10]. We follow 
the standard convention of putting uji = cDq, that is 

UJi = mf{u : ||T;|| < ikLe""* for all t > 0}. 

Theorem 33 For any < cr < 1 there exists a one-parameter semigroup Tt acting 
on a Hilbert space Ti such that s = Q, uoq = 1 and uoi = a. 

Acknowledgements I should like to thank C J K Batty, R Nagel and L N Tre- 
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